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We review recent theoretical progress in evaluating higher order QCD corrections to Higgs 
^ \ boson differential distributions at hadron-hadron colliders. 

^ : 

^ 1 Introduction 

The origin of the spontaneous symmetry breaking mechanism (ssbm) in particle physics is still 
j> \ unknown. In the Standard Model (SM) [1] a vacuum expectation value (vev) v = 246 GeV is 
given to components of a complex scalar doublet. Three of these fields generate mass terms for 
the and Z bosons. The remaining scalar is called the Higgs [2] and has not been observed. 
The present lower limit for rrtu from the LEP experiments [3] is 114 GeV/c 2 . In the SM the 
interaction vertices with scalar fields are expressed in terms of v and m#. Extensions of the 
standard model to incorporate supersymmetry contain several Higgs particles, both scalar and 
pseudoscalar [4]. In the minimal supersymmetric extension of the Standard Model (MSSM), there 
are two Higgs doublets, so it is called the Two-Higgs-Doublet Model (2HDM). After implementing 
the ssbm there are five Higgs bosons usually denoted by h, H, (CP even), A (CP odd) and H ± . 
At tree level their couplings and masses depend on two parameters, and the ratio of two 
vevs, parametrized by tan/? = v%/v\. Since no Higgs particle has been found the region < 
92 GeV/c 2 and 0.5 < tan/? < 2.4 is experimentally excluded [5]. For recent reviews on the 
theoretical status of total Higgs boson cross sections see [6]. Let us denote the neutral bosons 
h,H,A collectively by B. It is possible that a B will be detected at the Fermilab Tevatron p — p 
collider (\/S = 2 TeV). If not it should hopefully be found at the Large Hadron Collider (LHC), 
a p — p machine under construction at CERN (\/S = 14 TeV). If no B is found then the 
ssbm must be realized in a different way, possibly via dynamical interactions between the gauge 
bosons. In parallel with the ongoing experimental effort higher order quantum chromodynamic 
(QCD) corrections to B production differential distributions are needed. The leading order (LO) 
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partonic reactions were calculated quite a long time ago and the next-to-leading order (NLO) 
corrections to these distributions only recently. Note that the detection of the B via its decay 
products depends on m B . The present mass limits allow the decays B — > 77, B — > bb and 
B — > WW*, ZZ* where W* and Z* are virtual vector bosons, which are detected as leptons 
and/or hadrons (jets). As itlb increases other decay channels (for example B — > W + W~ , ZZ , 
B — > tt), open up requiring different experimental triggers. Higher order QCD corrections to 
these decay rates have also been calculated but will not be discussed here. 

2 Higgs differential distributions 

For inclusive B-production one calculates the differential cross section in the B transverse mo- 
mentum (jpx) and its rapidity (y), which are functions of m^, the partonic momentum fractions 
Xi, x 2 in the hadron beams, and the partonic cm energy y/s = ^xix 2 S. All other final state 
particles are integrated over. In contrast exclusive B-production retains the information on these 
other particles. In the SM the H couples to the gluons via quark loops with the qqH vertex pro- 
portional to m q , so the t-quark loop is the most important. In the 2HDM the ggA amplitude 
with quark loops depends on both the masses of the quarks and f3. In LO the g + g — > H cross 
section (order a 2 ) containing the top-quark triangle graph, was computed in [7]. However here 
Pt = 0, so we need a two-to-two body partonic process to produce a Higgs with a finite p T . 
Note that these are NLO processes with respect to the total B production cross section and of 
order a^. At small x the gluon density g(x, Q 2 ) > qi(x, Q 2 ) > q~i(x, Q 2 ), (?j stand for u, d, s, c, b, t 
quarks) so we expect that, in order of importance, the dominant production channels are 

9 + 9 ^9 + B, 
9 + ft(ft) -> ft(ft) + B, 

q i + q i ~^9 + B. (2.1) 

These are the LO order Born reactions for Higgs px and y distributions. The corresponding 
Feynman diagrams contain heavy quark box graphs. The B differential distributions for the 
reactions in Eq. (2.1) were computed for B=H in [8] and for B=A in [9]. The total cross section, 
which also contains the virtual QCD corrections to the ggB top-quark triangle, was calculated in 
[10], [11] and [12]. The expressions in [12] for the two-loop graphs with finite mt and tub are very 
complicated. Furthermore also the two-to-three body reactions (e.g. g + g—>g + g + B involving 
pentagon loops) have been computed in [13] using helicity methods. From these results it is 
clear that it will be very difficult to obtain the NLO (order af) corrections to the B differential 
distributions as functions of both m t and mg. 

Fortunately one can simplify the calculations if one takes the limit m t — > 00. In this case the 
Feynman graphs are obtained from an effective Lagrangian describing the direct Bgg coupling. 
An analysis in [14] in NLO reveals that the error introduced by taking the m t — > 00 limit is less 
than about 5% provided m B < 2 m t . The two-to-three body processes were computed with the 
effective Lagrangian approach for the B in [15] and [16] respectively using helicity methods. The 
one-loop corrections to the two-to-two body reactions above were computed for the H in [17] 
and the A in [9]. These NLO matrix elements (order af) were used to compute the pt and y 
distributions of the H in [18], [19], [20], [21] and the A in [22]. For ways to differentiate between 
production of H and A see [23] and references therein. 
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In the large m t limit the Feynman rules (see e.g. [15]) for scalar H production can be derived 
from the following effective Lagrangian density 



j0^ ff = G u ^ u (x)O(x) with 0{x) = --G%{x)G a ^{x), (2.2) 
whereas pseudoscalar A production is obtained from 

Ctf f = ^ a (x){G a 1 (x) + G A 2 (x)) , with 

i i ri f 

1 (x) = --e^G^(x)G^(x) , 2 (x) = --d»J2Vi(x)l»lsQi(x), (2-3) 
8 1 i=i 

where $ H (x) and $ A (x) represent the scalar and pseudoscalar fields respectively and rif denotes 
the number of light flavours. Furthermore the gluon field strength tensor is given by G^ u (x) (a 
is the colour index) and the quark field is denoted by <2i(x). The factors in the definitions of 
0,Oi and 2 are chosen in such a way that the vertices are normalised to the effective coupling 
constants G H) G A and G A . The latter are determined by the triangular loop graphs, which 
describe the decay processes B — > g + g with B = H, A, including all QCD corrections and taken 
in the limit m t — > oo, namely 



G, 



-2 5/4 a s (rf)G)( 2 T B F B (T B )C B \aM,£jj , 



G A 



where a s ([j%) is defined by 



a s (rf)C F (I -3 In 4 J + 



G A , 



(2.4) 



47T 



(2.5) 



with ck s (// 2 ) the running coupling constant and \i r the renormalization scale. Further Gp repre- 
sents the Fermi constant and the functions F B are given by 



F H (r) = 1 + (1 - r) f(r) , F A (r) = f(r) cot (3 , 
/(r) = arcsin 2 , for r > 1 , 



4 m 2 



m 



B 



/ M = -i( ln TTvS + " / 



for t < 1 , 



where (3 denotes the mixing angle in the 2HDM. In the large m t -limit we have 

lim F h (t) = — , lim F a (t) = - cot (3 . 

T— >oo O T r^oo T 



(2.6) 



(2.7) 



The coefficient functions C B originate from the corrections to the top-quark triangle graph pro- 
vided one takes the limit m t — > oo. The coefficient functions were computed up to order a 2 s in 
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[14], [24] for the H and in [25] for the A. The answer depends upon af\ the running coupling in 
the five-flavour number scheme, and we only need the first order terms 

C„(a s (/i r 2 ),4) = l + lla?> + ..., 

v mf ' 

C A (a s (fi 2 r ),^)=l. (2.8) 
The last result holds to all orders because of the Adler-Bardeen theorem [26] . 



3 Numerical results at moderate pr 

The hadronic cross section da for Hx(Px) + H 2 {P2) B(—P 5 ) + X is obtained from the partonic 
cross sections da a b for the reactions in Eq. (2.1) and their NLO corrections [for example g(pi) + 
9 (pi) -> g{-Pz) + g{-Pi) + B(-p 5 )] using 

(11 (1 U a,b=q,g Xl ' min 1 Jx 2,min -1-2 

x / b H2 (x 2 , /x 2 ) s 2 -^- (a, t, «, ml, /i 2 ) . (3.9) 

Here f^(x, fi 2 ) is the parton density for parton a in hadron H at factorization/renormalization 
scale \x and the hadronic kinematical variables are defined by 

S = (Pi + P 2 ) 2 , T=(P 1 + P 5 ) 2 , f/=(P 2 + P 5 ) 2 . (3.10) 



The latter two invariants can be expressed in terms of the p T and y variables by 



T = m B — v SJpt + coshy + v SJpj, + mg sinhy , 



[/ = m| - \fS\Jp 2 T + m| cosh y - \fS\Jp^ + m| sinh y . (3-11) 

In the case parton pi emerges from hadron Hi(Pi) and parton p 2 emerges from hadron if 2 (P 2 ) 
we can establish the following relations 

pi = x x Pi , £>2 = a; 2 P2 , 

s = iii2'5, i = %x(T — m 2 B ) + , u = x 2 (U — m 2 B ) + m\ , 



S • T m 2 ; • 2jni: ' .r,.S • T -m| • ' '"' 



Since the differential cross section contains terms in m B /p T In 1 (tub /pt), (i — 1? 2, 3) which are 
not integrable at pt = 0, we cannot integrate over p T down to pt = to find the y-distribution. 
However we can integrate over y to find the pr-distribution, valid over a range in pt where there 
are no large logarithms, (say pr > 30 GeV/c), 

d a HlHa ,„ o , p-ax , rf2 (J H 1 H 2 



(S,PT,m? B )= f dy ° (S,PT,y,m 2 B ) , (3.13) 
a J-?/max dp T dy 
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Figure 1: The differential cross section d a/dpx integrated over the whole rapidity range (see 
Eq. (3.13)) with m fl = 120 GeV/c 2 and fi 2 = m% + p\. The LO plots are presented for the 
subprocesses gg (long-dashed line), q(q)g (dot-dashed line) and 100 x (qq) (dotted line) using 
the parton density set MRST98(lo05a.dat). 

with a fixed y max . The calculation of the NLO B differential distributions requires the virtual 
corrections to the reactions in Eq. (2.1) and the NLO two-to-three body reactions (order oQ). 
Hence one needs a regularization scheme, and renormalization and mass factorization (which 
introduces the scales ji r and /!/ respectively). The presence of the 75 matrix in the pseudoscalar 
contributions makes everything even more complicated and we refer to [22] for details. In [18] 
helicity amplitudes were used and the fully exclusive two-to-three body reactions were calculated 
numerically. Very few formulae were presented. In [19] the cancellations of the UV and IR 
singularities were done algebraically leading to the H inclusive distributions. Many analytical 
results were given but the two-to-three body matrix elements were too long to publish. In [20], 
[21] helicity amplitudes were used for the H inclusive calculation and complete analytical results 
were provided. The numerical results from these three papers have been compared against each 
other and they agree. 

We put rif = 5 in a s (/i 2 ), <j a b and f^(x,fi 2 ) in Eq.(3.9). For simplicity /i r — /if — /i and we 
take /i 2 = m 2 H + p\ for our plots. Further we have used the parton density sets MRST98 [27], 
MRST99 [28], GRV98 [29], CTEQ4 [30] and CTEQ5 [31]. 

We want to emphasize that the magnitudes of the cross sections are extremely sensitive to 
the choice of the renormalization scale because the effective coupling constants in Eq.(2.4) are 
proportional to a s (/i r ), which implies that da LO ~ and dcx NLO ~ aj. However the slopes of 
the differential distributions are less sensitive to the scale choice if they are only plotted over a 
limited range. For the computation of the ggB effective coupling constants in Eq. (2.4) we take 
m t = 173.4 GeV/c 2 and Gp = 1.16639 GeV~ 2 = 4541.68 pb. Here we will only give results for 
H production at the LHC. Lack of space limits us to showing only p^-distributions. The LO and 
NLO pr differential cross sections in Eq. (3.13) with mn = 120 GeV/c 2 are shown in Figs. 1 
and 2 respectively. The MRST98 parton densities [27] were used for these plots. We note that 
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Figure 2: Same as Fig. 1 in NLO except for 100 * abs(gg) (dotted line) and the additional 
subprocess 100 x abs(gg) (short-dashed line) using the parton density set MRST98(ft08a.dat). 

the NLO results from the q(q)g and qq channels are negative at small pt so we have plotted 
their absolute values multiplied by 100. Clearly the gg reaction dominates. The q(q)g reaction 
is lower by a factor of about five. 

Regarding the corresponding distributions for A-production they have the same shape depen- 
dence in LO and differ slightly in NLO. Therefore the only difference is in the overall couplings 
in Eq. (2.4). In particular if tan/5 = 1 then da a/ dan = 9/4 for all values of rnn and itla in LO. 
There are small differences from 9/4 in NLO. Plots are presented in [22]. 

We now show the scale dependence of the distributions. We have chosen the scale factors 
\i = 2/i , \i = /i and /i — /iq/2 with /Xq = m 2 H + p\ and plot in Fig. 3 the quantity 



in the range 0.1 < /j/a*o < 10 at fixed values of p T = 30, 70 and 100 GeV/c. The upper set of 
curves at small /i//io are for LO and the lower set are for NLO. Notice that the NLO plots at 
70 and 100 are extremely close to each other and it is hard to distinguish between them. One 
sees that the slopes of the LO curves are larger that the slopes of the NLO curves. This is an 
indication that there is a small improvement in stability in NLO, which was expected. However 
there is no sign of a flattening or an optimum in either of these curves which implies that one 
will have to calculate the differential cross sections in NNLO to find a better stability under 
scale variations. 

Next we show the mass dependence of the NLO result in Fig. 4 using the MRST99 parton 
densities. The differential distribution drops by a factor of two as m H increases from 120 to 180 



There are two other uncertainties which affect the predictive power of the theoretical cross 
sections. The first one concerns the rate of convergence of the perturbation series which is 



N (pr, ± 



da(p T ,fi)/dp T 
da(p T ,/i )/dp T 



(3.14) 



GeV/c 2 . 
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Figure 3: The quantity N(p T: /-f/Vo) (see Eq. (3.14)), plotted in the range 0.1 < ///Vo < 10 
at fixed values of p T with m H = 120 GeV/c 2 and /!q = + using the MRST98 parton 
density sets. The results are shown for px = 30 GeV/c (solid line), pr = 70 GeV/c (dashed 
line), pt = 100 GeV/c (dot-dashed line). The upper three curves on the left hand side are the 
LO results whereas the lower three curves refer to NLO. 



indicated by the if-factor defined by 



K 



d a 



NLO 



da LO 



(3.15) 



Depending on the parton density set the i^-factors are pretty large and vary from 1.4 at pt = 
30 GeV/c to 1.7 at pt = 150 GeV/c for both H and A production. Another uncertainty is 
the dependence of the px distribution on the specific choice of parton densities, which can be 
expressed by the factors like 



^>CTEQ _ 



d a MRST 



j^GRV _ 



d a 



GRV 



d a MRST ' 



(3.16) 



and are generally above unity. Again these factors are essentially identical for H and A produc- 
tion. For specific values consult [18], [19], [21] and [22]. 

The reason why the parton density sets yield different results for the pr-distributions can be 
mainly attributed to the small x-behaviour of the gluon density because gluon-gluon fusion is the 
dominant production mechanism. Future HERA data will have to provide us with unique gluon 
densities before we can make more accurate predictions for the Higgs differential distributions. 



4 Numerical results at large x. 

Near threshold the longitudinal momentum fractions approach unity. In this region the soft- 
plus- virtual (S+V) gluons and collinear q — q pairs dominate the NLO corrections to the partonic 
cross sections. The S+V gluon parts of these cross sections are obtained by omitting the hard 
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Figure 4: The mass dependence of d <t nlo jdpx (see Eq. (3.13)) using the set MRST99(cor01.dat) 
with /i 2 = m 2 H +Pt for Higgs masses mn = 120 GeV/c 2 (solid line), m H = 160 GeV/c 2 (dashed 
line) and m H = 200 GeV/c 2 (dot-dashed line). 

contributions which are regular at s 4 = (p 3 + p 4 ) 2 = and adding the pieces from the virtual 
contributions. These two contributions constitute the S+V gluon approximation. To study its 
validity we show in Fig. 5 the ratio 

dS+V _ a a (a -1 7 \ 

1X J ^-EXACT ' V ' > 

for the NLO contributions to the px distribution (here px.min = 30 GeV/c). One expects that the 
approximation becomes better at larger transverse momenta where px approaches the boundary 
of phase space at x — 1. However in Fig. 5 the highest value of px, given by px = 150 GeV/c, is 
still very small with respect to PT,max ~ V / 5'/2 = 7 x 10 3 GeV/c. Therefore it is rather fortuitous 
that the approximation works so well for p T > 100 GeV/c where one obtains R s+Y < 1.2. 

The S+V gluon approximation overestimates the exact NLO result but the difference de- 
creases when the pt increases. In particular for px > 200 GeV/c the S+V approximation is 
good enough so that resummation techniques could be used to give a better estimate of the 
Higgs boson p T distribution corrected up to all orders in perturbation theory. Note that the 
boundary of phase space is also approached when m B increases at fixed px, but, according to 
the results in Fig. 5, the S+V approximation does not improve. In [21] it was noted that in- 
creasing tub makes the terms in Xn{mB/pT) larger at fixed px- So the range in px where the 
small ^-logarithms dominate increases with increasing thb- This leads us to the next topic. 
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Figure 5: The ratio i? s+v in Eq. (4.17) for the pt distributions using the set MRST99(cor01.dat) 
with jj? = m 2 H + PT,mm an d various Higgs masses given by mu = 120 GeV/c 2 (solid line), 
m H = 160 GeV/c 2 (dashed line) and m H = 200 GeV/c 2 (dot-dashed line). 

5 Numerical results at small pr 

We noted previously that there are terms in uib/pt ^(ttib/pt) which are dominant in the small 
Pt region. In this region the double differential cross section can be expanded as follows 

^ = 4[EE(?)"e..(4)* 

dp T dy s p T rrv^i ^ V 2vr^ V p 2 J 

where the next order terms start with m = 3 (order a^), and <r denotes one of the partonic 
cross sections for the LO processes in Eq.(2.1), which are order a 2 . In [21] the C mn were 
determined from their explicit NLO results in terms of certain factors A^ 2 \ B^ 2 \ 

C^J, CfV which multiply convolutions of splitting functions with parton densities. These factors 
had been previously determined in [32], [33] from other reactions, so it was gratifying to see 
the consistency between the different results. In Fig. 6 we show a comparison of the inclusive 
result in Eq. (3.9) at both LO and NLO versus the corresponding result from the small p^-limit 
formula in Eq. (5.18). The latter works very well for p? smaller than 10 GeV/c and agrees with 
the exact calculation for p? smaller than 5 GeV/c. The authors then propose another large x 
approximation involving both the S+V and small pt terms. We refer to their paper for details. 

Note that the resummation of the logarithms in Eq. (5.18) can be carried out following the 
procedure in [34]. This leads to a change in the shape of the p^-spectrum at small pt- Results 
can be found in [35], [36], [37]. 

In conclusion we note that the NLO corrections to the B differential distributions are now 
completely known and resummation methods have been used to study the region near p T = 0. 

Acknowledgement: We thank C. Glosser and C. Schmidt for giving us the inputs for Fig.6. 
We also thank B. Field, M. Tejeda-Yeomans, S. Dawson R. Kauffman, D. de Florian, M. Grazzini 
and Z. Kunszt for discussions and comments. 
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Figure 6: The Higgs p T spectrum compared to the small-p T limit formula (see Eq. (5.18)) at 
both LO and NLO with CTEQ5L and CTEQ5M1 parton densities respectively. All curves are 
calculated for y = in the m t = oo effective theory. 
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